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Abstract 

Finding all Bell inequalities for a given number of parties, measurement settings, and mea- 
surement outcomes is in general a computationally hard task. We show that all Bell inequalities 
which are symmetric under the exchange of parties can be found by examining a symmetrized 
polytope which is simpler than the full Bell polytope. As an illustration of our method, we gen- 
erate 238885 new Bell inequalities and 1085 new Svetlichny inequalities. We find, in particular, 
facet inequalities for Bell experiments involving two parties and two measurement settings that 
are not of the Collins-Gisin-Linden-Massar-Popescu type. 

1 Introduction 

Already discovered by Boole in the theory of logic and probabilities as "conditions of possible 
experience" [1], Bell inequalities found a new dimension with the work of John Bell who showed 
that quantum physics could violate these conditions in some situations, highlighting what is now 
known as quantum nonlocality [2]. 

Complete set of Bell inequalities are known only for setups involving small numbers of parties, 
measurement settings, and measurement outcomes. This may already be sufficient for various 
applications, such as exhibiting the nonlocality of a noisy quantum state in a real experiment [3], or 
establishing the security of a device-independent quantum key distribution protocol jUIS]. But the 
simplest inequalities are not always optimal. For instance, certain inequalities with large number 
of measurement settings are much more resistant to the detection inefficiencies than the CHSH 
inequality [U [7] , or are violated by quantum states that do not violate the CHSH inequality [U [9] . 
This motivates a search for Bell inequalities involving more parties, measurements, or outcomes. 
But finding all Bell inequalities pertaining to a given experimental setup is in general a hard task 
|XQ(, lllj . and a complete search with current-day techniques is not feasible in most instances. 

It is instructive, however, to realize that many useful Bell inequalities, like CHSH [12], Mermin 
|13j or CGLMP [14J to cite just a few, can be written in a form that is invariant under any 
permutation of the parties. Symmetric inequalities are also attractive because they are likely to be 
easier to handle. Motivated by this observation, we show here how to exploit a symmetric version 
of the full Bell polytope to generate all symmetric Bell inequalities. This symmetrized polytope is 
much easier to handle than the full Bell polytope. In particular, for Bell experiment with binary 
settings and outcomes, the number of extremal points of the symmetrized polytope only grows 
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polynomially with the number of parties, which is an exponential gain compared to the general 
situation. 

Our method for finding symmetric inequalities is not restricted to Bell inequalities, but can 
be applied to any set of inequalities characterizing a correlation polytope, like for instance the 
Svetlichny inequalities, which allow to test for genuinely multipartite nonlocality [15j. We present 
in the next section our approach from this general perspective. We then apply in Section 3 our 
method to several examples for which listing all (non-symmetric) inequalities is computationally 
intractable with present-day techniques. In the Bell scenario with two parties, two settings, and 
four outcomes, we find in particular facet inequalities that are not of the CGLMP form, answering 
an open question raised by Gill [T6] . 



2 General setting 

Let (n, m, k) denote a Bell experiment where n parties can choose one out of m possible measure- 
ment settings that each yield one out of k possible outcomes^} The statistics of the observed results 
are described by the joint conditional probability distributions (also called correlations^]) 

p(ri,...,r n \si,...,s n ) (1) 

where Sj £ {1, . . . , m} denotes the measurement setting of party i and n E {1, . . . , k} denotes the 
corresponding measurement result. 

Note that in general the N = m n k n probabilities ([I]) are not all independent but satisfy linear 
constraints, such as the normalization conditions or the no-signalling conditions. We are thus 
actually interested in an affine subspace of M> N of dimension d, where d = m n (k n — 1) for normalized 
correlations, and d = {m{k — 1) + l) n — 1 for correlations that satisfy in addition the no-signalling 
conditions [91 [T7] . We suppose in the following that a proper parametrization has been introduced 
so that the joint distributions ([!]) can be identified with points p in M. d . 

We are interested in whether a given p belongs to some special subset P C M. d . In this work, 
we consider sets P which are polytopea^J A polytope can be described by the list V of its vertices 
(or extremal points) v GV, and any point p 6 P can be written as 

P=^PvV, (2) 

where p v are positive and normalized weights: p v > for all v and Y^veV P« = 1- I 11 ^ ne case °^ * ne 
Bell polytope, for instance, the extremal points are the deterministic local strategies, corresponding 
to the prior assignment of an outcome T%& iV to each measurement setting Sj. There are thus k nm 
different vertices v, corresponding to joint probability distributions of the form 

/ i J 1 if rj = r i Si V for alH = 1, . . . , n 

p{n, . . .,r n \si, . . . ,s n ) = < (3) 

I U otherwise . 

1 More general Bell scenarios with different number of measurement settings mi and outcomes hi for each party i 
are also possible, but since we will consider situations that are symmetric under permutations of the parties, we 
choose m,i = m, ki = k for all i. 

2 Throughout this paper, the term "correlations" refers to probability distributions of the form It should not 
to be confused with "correlator" or "correlation functions", such as E(si,S2) = p(ri = r2|si,S2) — p[ri / r2|si,S2). 

3 In general, it might also be interesting to consider sets that are not polytopes |18U19| . nor even convex [20] , 
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Alternatively to its representation in term of vertices (the y-repesentation) , a polytope P can 
also be described, by the Farkas-Minkowski-Weyl theorem, as the intersection of finitely many half- 
spaces (the ^-representation) . A half-space is defined by an inequality h ■ p = Yli=i^iPi — 
specified by the couple (h, ho) G We say that an inequality (h, ho) is valid for the polytope 

P, if h- p < ho for all points p in P. The Farkas-Minkowski-Weyl theorem states that a polytope 
can be characterized by a finite set of valid inequalities. That is, a point p belongs to P if and 
only if 

h ■ p < h for all (h, h ) G H , (4) 

where H is some finite set in This description is particularly appropriate when willing to 

show that a point does not belong to the polytope, since it is sufficient to exhibit the violation of a 
single one of the inequalities Q. In the case of the Bell polytope, these inequalities are known as 
Bell inequalities. 

A complete and minimal representation of a polytope in the form Q is provided by the set 
of facets of the polytope. An inequality (/, fo) G H defines a facet if its associated hyper-plane 
/ • p = fo intersects the boundary of the polytope in a set of dimension d — 1, i.e., if there exists 
d affinely independent points of P satisfying f ■ p = fo- The set f = {p\p€.P 1 f-p = po} then 
corresponds to the facet defined by (/, fo). 

Finding all the Bell inequalities corresponding to an experimental configuration (n, m, k) thus 
amounts to determining the facets (minimal H-representation) of a polytope when given its extremal 
points (V- representation). This conversion problem is well known and there exists several available 
algorithms to solve it |21j. However, when n, m, or k are too large, the associated polytope becomes 
too complex to be handled by these algorithms. It might then be appropriate to focus the search 
on a subclass of all facet inequalities. We explain how this can be done for symmetric inequalities 
in the next section. 

2.1 Focusing on symmetric facets 

Some facet inequalities (/, fo) of a correlation polytope P can be written in a way that is invariant 
under permutations of the parties. This is the case for instance for the CHSH inequality, which 
can be written in the CH form as 

p(oi) + p(h) - p(aih) - p(aib 2 ) - p(a 2 h) + p{a 2 b 2 ) > (5) 

where we write p(a,b\x,y) for p(r±, r 2 \s±, s 2 jft and define p(a x ) = p{a = l\x), p(b y ) = p(b = l\y), 
p{a x b y ) = p{a = 1,6 = l\x,y). In the following we will call such facets symmetric facet^ If one 
is interested in finding only the symmetric Bell inequalities relevant to a given scenario (n,m,k), 
then it is possible to restrict the space in which to search for them. This is what we show now. 
Let G be the set of permutations of {1, . . . , n}. Given a permutation tt G G let 

p(ai, . . . , Or, | XI, . . . , X n ) l-> P(a7r(l), • • • «7r(n)l^7r(l), ■ • • ! ^vr(n)) • (6) 

be its action on the joint probability distributions. This permutation induces a transformation 
7r : M> d i— > ~R d : p i— > up in the vector space ~ML d in which the correlations p are represented, which 

4 From now on, we write p(a, b,c. . .\x,y, z, . . .) for p(ri, r2,rs, . . . \s%, S2, S3, . . .). 
Note that a symmetric inequality need not appear in a symmetric form when written in any of its equivalent 
forms under relabeling of measurement settings and outcomes. For instance : p(ai) + p(&2) — p(ai6i) — p(a\b2) + 
p(a2&i) — p(a2&2) > is equivalent to the CH inequality |5} up to relabeling of the settings and outcomes, but is not 
invariant under the exchange of party 1 and 2. 
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by abuse of language we denote by the same symbol ir. Note that the correlation polytopes that 
we consider here are evidently invariant under such permutations, i.e., irP = P, and any vertex 
v £ V is mapped into another vertex irv G V. 

Given a facet-inequality (/, fo) of P, let (irf, fo) be its image under ir. Note that with this 
definition the facet / = {p \ p G P, f • p = fo} is mapped onto wf = {irp \ p G P, f ■ p = fo} since 
{p I P G P, (tt/) • P = fo} = {P I P e P, f • (vrtp) = fo} = {up \ p G P, / • p = fo}, where we used the 
fact that the transformations ir are unitary, i.e., 7T = 7r _1 . We say that a facet (/, /o) is symmetric 
if 7r/ = / for all ir £ G. 

Consider now the symmetrizing map 

and let n denote the projection on the symmetric affine subspace S of M rf of dimension d s = dim(S'). 
We suppose to simplify the presentation that the origin of R d is contained in S (this can always be 
achieved by a proper translation of the correlation vectors p) , so that S is actually a linear subspace 
of M. d . An arbitrary vector p G M. d can then be written as p = p s +pt where p s = Yip is the projection 
of p on the symmetric subspace S and pt = (1 — n)p on the complementary space T = (1 — IT)IR d . 
Similarly an arbitrary inequality (/, /o) can be written as (f s © ft,fo)- A symmetric inequality 
then takes the form (f s © 0, fo) where denotes the null vector in ~R. d ~ ds . Given an inequality 
(fs,fo) £ M ds+1 defined in the symmetric subspace S, we denote its symmetric extension as the 
inequality (/, fo) = (f s © 0, fo) defined in the full space R d . 

Let P s = IIP = {p s | p G P} be the projection of the polytope P on the symmetric subspace. 
Any vertex w G P s is necessarily the projection w = v s of some vertex v G P. Indeed, suppose 
that w is the projection of a non-extremal point p = qiVi G -P, % = 1, < qi < 1. Then it 
needs to be the projection of every vi as well: w = lip = ^ qiflvi = w implies w = Tlvi since w 
is a vertex. Note, however, that every vertex v G P does not necessarily induce a vertex v s G P s 
when projected on the symmetric subspace. The symmetrized polytope P s has thus in general less 
vertices than the original polytope. Moreover it is defined in a space of smaller dimension d s < d 
than the full space M. d . It is thus in general easier to determine the facets of the symmetrized 
polytope P s than those of the full polytope P. The following theorem shows that determining the 
facets of this symmetrized polytope is sufficient to find all symmetric facet inequalities of the full 
polytope P (see also Figure 1). 

Theorem. Let (f s , fo) be a facet inequality for the polytope P s G ~R ds . Then its symmetric extension 
{fs © 0, fo) G M. d+1 defines a valid inequality for the full polytope P G Moreover all symmetric 
facet inequalities of the full polytope P are the symmetric extension of some facet of the symmetrized 
polytope P s . 

Proof. The symmetric extension (/, fo) = (/ s © 0, fo) is valid for P if / ■ p < fo is satisfied by all 
points p G P. But this immediately follows from the fact that / -p = (f s ©0) ■ (p s ®pt) = fs'Ps < fo 
and the fact that f s -p s < fo is valid for all p s G P s . 

Now, let (g,go) = (9s © 0,go) be a symmetric facet of P. Clearly, it is the symmetric extension 
of the inequality (g s ,go) G M. ds , which is valid for P s . Moreover, (g s ,go) defines a facet of P s as 
there exist d s affinely independent points in P s that saturate it. Indeed, since (g, go) defines a facet 
of P, there exists d affinely independent points p in P that satisfy (g s © 0) • p = go- These points 
are of the form p = p s ® pt, where pt can clearly be arbitrary. Since the complementary space T is 
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Figure 1: (a) Example of a polytope P in the vector space M 3 . (b) Subspace S symmetric under the 
exchange of coordinates e± and e 2 - P s (grey) is the projection of the polytope onto this subspace. 
(c) f s and g s are two facets of P s , and / and g are their symmetric extensions to the whole space 
M 3 . / is a symmetric facet of the original polytope P, whereas g is just a valid inequality for P. 



of dimension d — d s , there must therefore be at least d — (d — d s ) = d s affinely independent points 
of the form p s © that saturate the inequality (g s © 0, go). These points obviously define d s affinely 
independent points in P s that saturate the inequality (g s ,go)- D 

Note however that the converse of the theorem is not true, as illustrated in Figure [TJ facets 
of P s do not necessarily extend to facets of the polytope P in the general space M d . We show in 
Section 2.3 how it is nevertheless possible to take advantage of such inequalities to generate new 
(not necessarily symmetric) facet inequalities for the original polytope P. 

2.2 Illustration on the (2, 2, 2) Bell scenario 

We now illustrate in detail the above approach on the (2,2,2) Bell scenario. This scenario is 
characterized by 16 probabilities p(ab\xy), where x S {1,2} denote the measurement setting of 
Alice and a G {1,2} the corresponding outcome, and where similarly y £ {1,2} and b £ {1,2} 
denote Bob's measurement setting and outcome. These probabilities satisfy normalization 

p(ab\xy) = 1 for all x, y = 1, 2, (8) 

a,6=l,2 

and no-signalling 

p(a\x) = p( a b\ x u) f° r an o, a;,j/ = 1,2 

6=1,2 

p(b\y) = p(ab\xy) for all b, x, y = 1, 2 . (9) 

a=l,2 

In total, only 8 of the 16 probabilities p(ab\xy) are therefore independent and we can represent the 
correlations p as elements of M 8 . For specificity, we choose the following parametrization 

P= [P(ai),p(a2),p(bi),p(b 2 ),p{aibi),p(aib2),p(a2b 1 ),p(a 2 b2)} , (10) 
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where p(a x ) = p(l\x), p(b y ) = p{l\y) and p(a x b y ) = p(ll\xy). 
The Bell-local polytope is described by 16 vertices 

Vaia 2 bi& 2 = [ai,a2,bi,b2,aib 1 ,a 1 b 2 ,a2bi,a2b 2 ] 



(11) 



where 01,02,61,62 € {1,2} specifies the deterministic assignment of an outcome to each measure- 
ment setting. 

The group G of permutation of two parties contains 2 elements: the identity 1 and the permu- 
tation 7r acting as follows on a vector p, 



7T p= [p(b 1 ),p(b 2 ),p(ai),p(a2),p(a 1 bi),p(a2bi),p(a 1 b 1 ),p(a 2 b2)] ■ 



(12) 



The symmetrizing map projecting on the space S of symmetric correlations is thus equal to II = 
|(1 + 7r), while the map projecting on the complementary space T is 1 — II = i(l — ?r). Arbitrary 
correlations p can thus be decomposed into a symmetric and an asymmetric part p = p s ®Pt, where 

p + np 



P, 



p(ai)+p(h) p(a 2 )+p(6 2 ) p{a 1 )+p(b 1 ) p(a 2 ) + p(b 2 ) 



p(aih), 



2 2 2 2 

p(aib 2 ) + p(a 2 bi) p(aib 2 ) + f>(a 2 6i 



(13) 



— ,p(a 2 b 2 



and 



Pt 



p — TTp 



p(ai)-p(bi) p(a 2 )-p(b 2 ) -p(ai) + p(h) -p(a 2 ) + p(b 2 



(14) 



p(aib 2 ) -p{a 2 bi) -p(aib 2 ) + p(a2h) 



Note that the symmetric part p s lives in a 5-dimensional subspace of M 8 and can thus be expressed 
in an appropriate basis as 



Ps 



p(ai) + p(b{) p(a 2 ) + p{b 2 ) 



,p(aibi 



p(aib 2 ) +p(a 2 b 



— ,p(a 2 b 2 



2 7 2 -v--/' 2 

Similarly, y>f lives in a 3-dimensional space of M 8 and can be decomposed in a proper basis as 

~p(ai)-p(bi) p(a 2 )-p(b 2 ) p(aib 2 ) - p{a 2 b 1 

Pt ~ 



(15) 



(16) 



The projection of the 16 deterministic points (11) on the symmetric subspace defined by (15) 
are given by 



Vs;a 1 a 2 b 1 b2 



a\ + 61 a 2 + 6 2 , ai6 2 + a 2 6i , 
ai6i , , a 2 6 2 



(17) 



2 ' 2 ' ' 2 
Note that some vertices of the original polytope are projected onto the same point of the symmetric 
polytope. For instance, v s -m2 = v S ;VZix- In total, it can be verified that the set defined by (17) 
contains only 10 extremal points. 



6 



We thus have reduced the original 8-dimensional polytope defined by 16 vertices to a 5- 
dimensional polytope with 10 vertices. Applying to this symmetrized polytope a standard al- 
gorithm performing the transformation from the V-representation to the H-representation, we find 
4 different types (up to relabeling of settings and outcomes) of facets of the symmetric polytope: 

P(aih) > 

P( ai b 2 ) + P(a 2 b 1 ) > 

P(ai) + P(h) - 2P(oi6i) > (18) 
P(oi) + P(h) - P(ai6i) - P(aife 2 ) - P(o 2 6i) + P(a 2 6 2 ) > 

We recognize the first inequality as the positivity condition for the joint probabilities and the 
last one as the CHSH inequality (written in the CH form as in equation ([5])). These two classes 
of inequalities define symmetric facets of the full polytope. The two other inequalities are valid 
inequalities for the full local polytope, but do not correspond to facets (although they are facets of 
the symmetrized polytope, as the inequality g in Figure [TJ. Note that in this (2,2,2) Bell scenario, 
the two only types of facet inequalities of the full polytope (the positivity condition and the CHSH 
inequality) can be written in a symmetric way. Hence in this simple case finding the facets of the 
symmetrized polytope is sufficient to generate all Bell inequalities. 

2.3 Generating facet inequalities from valid inequalities 

As we mentioned earlier, and was illustrated above, facets of the symmetrized polytope P s can 
correspond to inequalities which are not facets of the original polytope P. These inequalities are 
nonetheless valid inequalities which are satisfied by all points in P and which might be violated 
by points that do not violate any of the symmetric facets of P. These inequalities may be used to 
generate new (non-symmetric) facets of P. 

There exist various deterministic or heuristic algorithms which may generate new facet inequal- 
ities starting from a valid (not facet-defining) inequality, see for instance [25) . Here, we describe a 
procedure that can be used whenever the starting valid inequality corresponds to a high-dimensional 
face of P, i.e., when the number of affinely independent vertices that saturate the inequality is large. 
In this case, it is possible to find all the facets that contain this high-dimensional face by completing 
the list of vertices with all possible combination of vertices that do not saturate the inequality, as 
detailed by the following algorithm: 

1. Let (/, /o) be a valid inequality for P and let W = {v \ f ■ v = /o} be the set of vertices 
saturating this inequality. 

2. Let dim(W) denote the number of affinely independent points in W. 

• If dim(VF) = d, then (/, /o) is a facet of P. 

• If dim(W0 < d, let U = {v\ dim(W Uw)> dim(VT)}. 

For every u £ U, let (<?, go) be the hyperplane passing through the points in U. If (g, go) 
is a valid inequality for P, it now defines a face of P of dimension dim(VF) + 1; in this 
case, go back to point 1 with (g,go) as a starting inequality. 
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3 Applications 



We now illustrate our method in several situations for which generating the complete set of facet 
inequalities using standard polytope software [21] is too time-consuming to be feasible. Due to 
the large number of inequalities that we have found, we only explicitly write a few of them here. 
Complete lists of all the facet inequalities that we generated are posted on the website [22]. Our 
results are summarized in Table I. 

Note that we list here only inequalities that belong to different equivalence classes, where two 
inequalities are considered equivalent if they are related by a relabeling of parties, measurement 
settings, or measurement outputs, or if they correspond to two different liftings of the same lower- 
dimensional inequality [T7]. In appendix A, we introduce a parametrization of the correlation 
space that naturally induces several invariants for each equivalence class. These invariants are 
easily computed and are useful to determine quickly wether two inequalities are equivalent (two 
equivalent inequalities have equal invariants). 

3-1 (2,2,4) 

We first consider bipartite experiments with two measurement settings per sites and 4 possible 
outcomes. Note that the case (2, 2, 3) was completely solved by Kaszlikowski et al. [23] and Collins 
et al. [H], who showed that all facets of the (2,2,3) polytope either correspond to the positivity 
of probabilities, the CHSH inequality, or the CGLMP inequality p3]. The CGLMP inequality was 
introduced for any number of outcomes k > 3 in |14j . and Gill raised the question [16] whether all 
non-trivial facet inequalities of (2, 2, k) are of this form. 

Using our method, we found that the Bell polytope corresponding to (2,2,4) contains 12 in- 
equivalent symmetric Bell inequalities. Among them, 8 involve the four possible outcomes in a 
nontrivial way, i.e., they correspond to genuine 4-outcome inequalities that cannot be seen as lift- 
ings of inequalities with lower numbers of outcomes, these inequalities are listed in Appendix B. 
The list of these 8 inequalities contains the CGLMP inequality, but surprisingly, it also contains 7 
inequalities that are inequivalent to it, thus answering in the negative Gill's question [16J. 



Bell scenario 


d 


d s 


\v\ 


\V S \ 


# symmetric inequalities 


# symmetric facets 


(2,2,4) 


48 


27 


256 


136 


29 


12 


(2,4,2) 


24 


14 


256 


136 


90 


55 


(4,2,2) 


80 


14 


256 


35 


627 


392 


(5,2,2) 


242 


20 


1024 


56 


>238464 


238464 


Correlators (3,3,2) 


27 


10 


512 


40 


44 


20 


Svetlichny (3,2,2) 


56 


14 


2944 


132 


1204 


1087 



Table 1: Summary of our numerical results. For each scenario we give the dimension of the space 
the polytope lives in as well as its number of extremal vertices, both before (d, \V\) and after 
projection on the symmetric subspace (d s ,|V^|). We give the number of inequivalent symmetric 
inequalities valid for P (but not necessarily facet-defining) obtained by resolving the symmetric 
polytope P s and the number of those inequalities that are facet defining. 
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3.2 (2,4,2) 



We now consider a bipartite scenario involving 4 settings with binary outcomes. The simpler case 
(2,3,2) was solved in |12 tl24l [9"] and contains a single new inequality besides the positivity constraints 
and the CHSH inequality. 

With 4 settings, we could use our method to find all of the 90 inequivalent facets of the sym- 
metrized polytope. Among these, 55 turn out to be facets of the (2,4,2) full local polytope, there 
are thus in total 55 symmetric inequalities for this scenario. Most of them were already known 
(see [TJ l25l 126]). but we could not find the two following ones in the litterature, given here in the 
notation of [9] : 



q51 

(2,4,2) 





-1 


-2 


-2 


-2 







-2 


-2 


-2 


-1 


-3 


3 


2 


2 





-3 


2 


-2 


1 


-2 


3 


2 


-1 


-1 < , 


c-52 _ „ 
J (2,4,2) — — * 


2 





2 


2 


-2 


2 


-1 


-1 


3 


-2 


-2 


2 


4 


-1 


-2 


2 


-1 


3 





-2 


1 


2 


-1 


1 



< 



(19) 



3.3 (4,2,2) and (5,2,2) 

Since our method takes advantage of the symmetry between parties, we expect that it will be 
particularly useful for multipartite Bell scenarios. Indeed for the Bell scenario (n,2,2), corresponding 
to n parties with binary settings and outcomes, the full local polytope has 4™ vertices and is 
embedded in a space of dimension 3 n — 1. The symmetrized polytope, on the other hand, has at 
most g(n+l)(n+2)(n+3) vertices and is embedded in a space of dimension d s = |n(n + 3). These 
quantities are polynomial in n and represent an exponential advantage with respect to the general, 
non-symmetric situation. Note that it therefore follows that it is possible using linear programming 
to decide in polynomial time in n if a given symmetric correlation vector p is local. 

The case (3,2,2) was already completely solved in [27]. For (4,2,2) we found a total of 627 
inequivalent symmetric inequalities, of which 392 are facet-defining. These facet inequalities corre- 
spond to the positivity conditions and to 391 genuinely 4-partite inequalities. Amongst them, the 
following one is quite interesting, as it can be violated by a 4-partite W state with measurements 
lying in the x-y plane (contrary to the 3-partite case where no inequality is known that can be 
violated by a W-state with measurements lying in the x-y plane): 



Iw = ~ p(aih) + p(aibici) + p(eti&ic 2 ) - p(a 2 b 2 C2) - p(ai&icidi) 

-p(ai&icid 2 ) - p(aihc 2 d 2 ) + p(aib 2 c 2 d 2 ) + p(a 2 b 2 c 2 d 2 ) + sym < 0. 



(20) 



The notation "sym" stands for the symmetric terms that are missing in (20 ), such as p(aiCi),p(oi6 2 ci) 
etc. If we consider the W state 1 0001) + 1 0010) + 1 0100) + 1 1000) and measurements in the x — y 
plane at an angle <p with respect to the x axis, and set 



<Pa 1 = 4>d 

<i>A 2 = 4>c 2 



o 



arccos 



>B 2 



1 



?D 2 



arccos 



2 arcsin 



-2 arcsin ■ 



1 



(21) 



we find a value Iw = 1/16 > 0. 
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For (5,2,2), we found 238464 inequivalent symmetric facets. Note that all of these inequalities, 
except the positivity of the probabilities, are truly 5-partite ones (i.e., they do not correspond to 
lifting of inequalities involving less parties). Among these inequalities, 9 of them involve only full 
(5-partite) correlators and were already given in [TT] . 

3.4 Correlation inequalities for (3,3,2) 

We considered also a tripartite scenario with 3 binary measurements per party. Since in this case 
even the symmetrized polytope is quite time-consuming to solve, we made a further restriction by 
considering only "full-correlator" inequalities, which can be written using only terms of the form 
(A x B y C z ) = p(a + b + c = 0|x, y, z) — p(a + b + c = l\x, y, z). This corresponds to performing a 
projection of the polytope on the subspace defined by 

{Ai) = (Bi) = {d) = (AiBj) = {AiCj) = {BiC 3 ) = for all i,j = 1, 2, 3. (22) 

We obtained 40 inequalities in this way, 18 of which are facets of the full original polytope that truly 
involve 3 inputs per party. Using the method presented in section [273| we found 13 supplementary 
facet inequalities, all of which involve again full-correlators only. 

3.5 Svetlichny inequalities for (3,2,2) 

To illustrate that our method is not restricted to Bell-local polytope, but can address any correlation 
polytope, we consider the Svetlichny polytope for 3 parties [15], which characterizes true tripartite 
nonlocality [281 l29l 130] . 

In a Svetlichny model, two of the three subsystems are allowed to communicate once the mea- 
surement settings have been chosen. There are thus three types of Svetlichny vertices uab/C) v ac/b-> 
v bc/Ai depending on which pairs of parties are allowed to communicate. A vertex of the form, e.g. 
v ABjCi corresponding to a deterministic strategy where outcomes a(x,y) and (3(x,y) are jointly 
determined for party 1 and 2, and an outcome j(z) is assigned to party 3. This defines a joint 
distribution of the form 

/ 7 I x J 1 if a = a(x,y),b = f3(x,y),c = 7(2) 
p{a,b,c\x,y,z) = < . (23) 

I otherwise . 

Such probability points do not satisfy the no-signalling conditions. For binary settings and out- 
comes, the Svetlichny polytope thus lives in a vector space of dimension d = 56. The subspace 
which is symmetric under the exchange of the three parties, however, has only dimension d s = 14. 

This great reduction in the space dimension, together with a reduction in the number of extremal 
points (see Table 1), allowed us to find all symmetric Svetlichny inequalities. After projection on 
the no-signaling space^ja total of 1087 facet symmetric Svetlichny inequalities were found. 

Interestingly, there are only two symmetric Svetlichny inequalities that involve only full-triparite 
correlation terms: the original Svetlichny inequality [TB] and the following one: 

Icorr = - {AxBxd) + {A 1 B 1 C 2 ) + {A x B 2 d) - 3(A 1 B 2 C 2 ) 

+ {A 2 B x Ci) - 3(^ 2j BiC 2 ) - 3(A 2 B 2 d) - 3(A 2 B 2 C 2 ) < 10 ( ' 

6 Indeed, we are only interested in whether these inequalities are violated by quantum correlations, which satisfy 
the no-signalling conditions. 
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This last inequality can be violated by quantum states, for instance by GHZ states having a visibility 
larger than 95.68%. 

The following inequality is also interesting: 

Ighz = -3P(a 2 ) + P(ai6 2 ) - P{a x bx Cl ) - P( ai b 2 c 2 ) + 7P(a 2 b 2 c 2 ) + sym < 0. (25) 

where "sym" stands for the missing symmetric terms. It can be shown [31] that it is violated by 
every GHZ-like state of the form 

\GHZ) = cos0|OOO) + sm0|lll). (26) 



4 Outlook 

Motivated by the number of interesting Bell inequalities that are invariant under permutations of 
the parties, we introduced a method to list all symmetric inequalities. This method works even 
in cases where solving the full correlation polytope is computationally intractable. Our method 
can also be used as a starting point to generate more general, non-symmetric inequalities using 



algorithms such as the one described in section 2.3 or in [25J. 



Our method allowed us to find a number of new Bell inequalities. But a new problem is now 
at sight: so many different inequalities are generated, even for simple situations, that it is difficult 
to find which ones are the most interesting. Evidence of this problem was already put forward in 
|llj and in [26] where it was shown that the number of inequivalent Bell inequalities increases very 
quickly with the number of parties or measurement settings. New insights are thus necessary in 
order to classify these inequalities and understand which ones are the most relevant. 

For simple Bell scenarios, such as (2,2,2), (2,2,3) or (2,3,2) all facet inequalities happen to 
be symmetric inequalities. What is the proportion of symmetric inequalities in more complicated 
scenarios? Are there other useful symmetries or properties that can be exploited to generate more 
inequalities? 
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Appendix A: Correlators and inequality invariants 



When dealing with Bell inequalities, it is useful to check quickly if two inequalities are equiva- 
lent under relabeling of parties, measurement settings, or outcomes. We were confronted with 
this problem when classifying the inequalities that we derived here. In this appendix, we intro- 
duce a parametrization of the correlation space which naturally induces several invariants for each 
equivalence class. Two inequalities that are equivalent have equal invariants. 

Let us consider a (n,m,k) Bell scenario and let I C {1, . . . ,n} be a subset of the n parties. 
Define p(ai\x) as the probability that these |/| parties obtain outputs aj = o^, . . . , a^, given that 
measurements x = x\, . . . ,x n have been made on all parties. Note that with this notation, the 
no-signalling condition is expressed as p(aj\x) = p(ai\xj), where xi = x^, ... , x^,. 

Define now single-party "correlators" E(ai\x) by 



E(a,i\x) = kp(ai\x) — 1 , 
and define by induction multipartite correlators E(aj\x) through 

E(ai 1 , ai 2 \x) = E{ai l \x) * E{ai 2 \x) 



(27) 



(28) 



where the * operation is just the usual multiplication, except when acting on two probabilities, 
in which case it satisfies p(aj 1 \x) * p(aj 2 \x) = p{aj 11 ai 2 \x). The no-signalling condition is then 
expressed as 

E(a T \x) = E(a I \x I ) . (29) 
The normalization condition on the probabilities p(ai\x) imply, on the other hand, that for any 

i e I 



£ 

Qi=0 



fc-1 

E(ai\x) = E(ai\x) * E(aj\i\x 

ai=0 
fc-1 

(kp(ai\x) - 1) 

cii=0 



(30) 



* E(a I \ i \x) = . 



The correlators E(aj\x) are in one-to-one correspondence with the probabilities p(aj\x) and 
thus represent an alternative parametrization of the correlation space. Note that in the case of 
binary outcomes {k = 2), E(ai\xj) coincides with the usual definition of a correlation function. 



The definitions (27) and (28) thus represent a possible generalisation of correlation function to 
more outcomes. 

With the notation that we just introduced, a generic Bell inequality in the no-signalling space 
takes the form 

c{a I: x I )E{a I \x I ) < c(0) (31) 

I,Xj,aj 

where c(o/, xj) are the coefficients of the inequality. A property of the correlators E(aj\xi) is that 
the white noise yields E{ai\xi) = 0. The resistance to noise of a an inequality written in the form 



(31) is thus directly given by the ratio of the local bound to the violation. 

A relabeling of parties, measurement settings, or measurement outcomes simply amounts to 
rearrange the order of the coefficients of the inequality. Note, however, that because of the nor- 



malization conditions (30), the basis of the correlation space that we chose is overcomplete. The 
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coefficients c(a/,x/) are thus not uniquely defined: adding J2 ai ^( a i\ x i) = f° r some i G I to the 
inequality ( |31[ ) does not change the inequality itself, but does change its coefficients. To compare 
two inequalities, we must therefore ensure first that they are written in some standard way. 



The freedom that we have in adding terms of the form A(i, an^xi) ^2 a . E(ai\xj) = to (31) 
corresponds to define new coefficients for the inequality in the following manner 



c'(aj, x£) = c(a/, xi) + X(i, atj^, zj) . 

We show now that requiring the inequality coefficients c'(ai,xj) to satisfy the relation 

fe-i 

c'(af , x/) = for all i £ J. 

<H=0 

allows to define them uniquely. 



(32) 



(33) 



Proposition. There exist values of X(i, an i5 xi) such that the newly defined coefficients c'(a/,x/) 
satisfy relation (33). Moreover for all such X's, the c'(ai,Xf) are the same, they are thus unique. 



Proof. Since equations (32) and (33) apply independently on every subset / of the parties, and on 



every inputs xj, we omit these indices, writing for instance c(ctj) instead of c(aj,xi) to lighten the 
notation. Moreover, all sums on the outputs ai go from to k — 1 and all sums on the parties i 
run on I, so we also omit these bounds in the proof. 

The existence of the A's can be shown by directly checking that the following formula is of the 



form (32), and satisfies (33). Let 



c'(ai) 



c{aj) 



(34) 



where we used the notation / + Ylaf + / = [1 + Sa + Sfe]/ f° r an y function / and the o 
composition satisfies 1 o 1 = 1, 1 o J2 a = J2 a ol = Y^ a , J2 a °^b = Safe an< ^ distributes over 
addition. To get this expression from (32) one possible choice of lambdas is: 



A(l,a Al ) 
A(2,a/\ 2 ) 



-r 



c(a/J 



a 1 



0,2 ai,a,2 



(35) 



and equation (33) is satisfied: 

E <='(»') = 



ai 



HE-E • 

\ 0>i 0>i / 



c(a/) = 0. 



c(a/) 



(36) 
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Now to show the unicity of the d coefficients, we notice that equation (33) is a non-homogeneous 
linear system of equations in the A variables, which we can write as: 



Y a A«) = - £ c ( a/ ) 



(37) 



Every solution of this system can thus be written as X = X p + X v where A p is a particular solution 
of the equation (as given by equation (35) for instance) and A^ is a solution of the homogeneous 
system, where the right-hand side of equation (37) is replaced by zero. Thus every d(ai) that 
satisfies equation ( 33 ) can be written as 

c'(ai) = c(ai) + ^ X p (i, a 7 \j) + ^ X v (i, aj\;). (38) 



Now we show that the last term of equation (38) is zero for every solution A^ of the homogeneous 



consider the following expression: 



counterpart of system (37), which implies that the coefficients d are uniquely defined. For this 



Z 



K(i,a>i\i)- 



It is clearly zero, since it contains the term 



X v (i,aj\i) = A v (i,a/\j) - ykX v (i, a^) = 0. 



On the other hand, we have that 



0.1 / 

Be 



o . . . o 



o . . . o 1 



"n-1 , 



O 1 



'-IS 

h 



^A„(i 5 o/\j) 

i 

^A v («,a/\j) 
^2X v (i,aj\i) 



l I\i) 



(39) 



(40) 



(41) 



(42) 

(43) 
(44) 



where the second term in (42) vanishes by definition of A„(i, ajw). Repeating iteratively the above 
step, we find eventually that 

= Y J Ui,a I \ i ). (45) 



This, combined with the fact that Z = 0, implies the desired result. 



□ 
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We showed that the coefficients of an inequality can be defined in a standard and unique way 
by requiring them to satisfy the constraints (33). Now, since a relabeling of parties, measurement 
settings, or outcomes can only rearrange the coefficients c'(ai,xi) without changing their value, the 
ordered lists of coefficients c'(aj, xi) for \I\ = 0, 1, . . . , n provide n+1 invariants for each equivalence 
class. For instance the local bound c(0) is an (easily checkable) invariant. In general requiring that 
two inequalities have their ordered list of coefficients identical does not guarantee that they are 
equivalent, but in the symmetric (4,2,2) scenario for instance, all the 391 classes of inequalities 
that we generated had different lists. 



Appendix B: List of symmetric inequalities with two inputs and 
four outcomes. 

The following inequalities are given in the notation of |9j. 

Case with 3 outcomes for the first setting and 4 outcomes for the second one. 



c2 

D (2,2,(3,4)) 





-1 


-1 











-1 





1 





1 


1 


-1 


1 





1 





1 








1 





-1 


-1 





1 





-1 





-1 





1 


1 


-1 


-1 


-1 



< 



(46) 
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Case with 4 outcomes for both settings. 



ol 

^(2,2,4) 



o3 

D (2,2,4) 



^(2,2,4) 



c-7 

^(2,2,4) 





-1 


-1 


-1 











-1 








1 


1 


1 


1 


-1 





1 


1 


1 


1 





-1 


1 


1 


1 


1 











1 


1 


1 


-1 


-1 


-1 





1 


1 





-1 


-1 








1 








-1 










-1 


-1 





-1 


-1 





-1 








1 


1 


1 





-1 





1 





1 





1 





1 








-1 


-1 





-1 


1 


1 


-1 


1 


2 





-1 


1 





-1 


2 


2 











1 











-1 




-1 


-1 


-1 











-1 





1 


1 








1 


-1 


1 





1 





1 





-1 


1 


1 





1 

















1 





-1 


-1 








1 





-1 





-1 





1 








-1 


-1 







-1 


-1 


-1 











-1 





1 


1 








1 


-1 


1 








1 


1 





-1 


1 





1 


1 














1 


1 


-1 





-1 








1 











-1 





1 








-1 


-1 






<o, 



< o, 



o2 

^(2,2,4) 



q4 

(2,2,4) 



<o, sf 



(2,2,4) 



<o, s: 



(2,2,4) 





-1 


-1 





-1 








-1 





1 





1 





1 


-1 


1 








1 


1 














-1 


1 








-1 


1 


1 


1 


1 














1 











-1 





1 











-1 







-1 


-1 





-1 


-1 





-1 


-1 


1 


-1 


1 


2 


1 


-1 


1 


1 





1 





1 





-1 





-1 


1 


1 





-1 


1 


1 


1 


1 








-1 


2 





1 





1 


-1 





1 


1 








-1 


-1 




-1 


-1 


-1 


-1 








-1 





1 


1 





1 





-1 


1 


1 





1 





1 


-1 


1 





-1 


2 





1 


-1 





1 


2 


1 


-1 








1 








-1 





-1 








1 


1 





-1 


-1 




-2 


-1 


-1 











-2 


2 





1 





1 


2 


-1 





1 


1 


1 








-1 


1 


1 








1 











1 





-1 


-1 








1 





1 


-1 





-1 





2 











-1 


-2 



< 



< 



(47) 
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